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Abstract. A self-consistent mathematical description of the 
magnetic field of fluted sunspot penumbrae is presented. This 
description is based on an expansion of the nonlinear force- 
free magnetohydrostatic equations written in cylindrical coor- 
dinates. The lowest order solutions are mathematically equiv- 
alent to laminated force-free equilibria in Cartesian geometry. 
The lowest order solutions have no toroidal component of the 
magnetic field and the magnetic pressure does not vary with 
azimuth but the solutions allow arbitrary variations of the mag- 
netic field components with azimuth. Explicit solutions are pre- 
sented which have a realistic radial profile of the magnetic field 
strength and reproduce the basic features of the observations. 
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1. Introduction 

Sunspot penumbrae show structures in the form of alternating 
dark and bright fibrils extending in the radial direction across 
the penumbra. Observations by Beckers & Schroter (1969) 
have already suggested that the magnetic field inside the dark 
filaments might be close to horizontal whereas the field in the 
bright filaments is less inclined with respect to the vertical. This 
magnetic field structure has implications for the Evershed flow 
which was also suggested to be associated mainly with the dark 
filaments. The association of the Evershed flow with the dark 
penumbral fibrils has recently been confirmed by Shine et al. 



(1994) and Rimmele (1994, 1995b). In these observations, it 



was also found that the Evershed effect is a quasi-periodic time- 
dependent phenomen on whi ch extends beyond the penumbral 
boundary (Rimm ele 1995b ). Furthermore, there is evidence 
(Rimmele 1995a) that the Evershed flow occurs in thin chan- 



nels which are elevated above the continuum height for most 
parts o f the penumbra. Very recently, Westendorp-Plaza et al. 
( 1997 ) have actually for the first time identified the regions of 
return flux for the field lines carrying the Evershed flow. Their 



observations support the view that the Evershed flow is carried 
by low-lying field lines. 

The suggested spatial variation of the inclination angle of 
the magnetic field with azimuth was also found by Lites et al. 



(199C) and Degenhardt & Wiehr (1991). Recently, very high 



resolution observations with the Lockheed tunable filtergraph 
on the Swedish Solar Observatory in La Palma (Title et al. 



1993) showed that the magnetic field of a sunspot penumbra 



indeed varies between nearly horizontal in the dark penumbral 
filaments and much less inclined in th e brigh t filaments. Similar 
results wher e obtained by Rimmele ( 1995a| ) and St anchfieldll 
et al. ( |1997| ). The observations of Title et al. ( |1993| ) show that 
the mean inclination of the magnetic field increases from 40°- 
50° to 70°-75° across the penumbra and that there is a rapid 
azimuthal variation of the inclination angle of about 18°. Little 
or no variation of the total field stren gth wi th azimuth is found 
by Title et al. ( 1993 ) and Rimmele ( 1995a ) (see, however, the 
discussion in Stanchfield ll et a l. ( 1997 )). On the basis of these 
observations, Title et al. ( 1993 ) proposed a tentative model for 
the magnetic field struc ture, which res embles that of laminated 
force-free fields (Low 1988a , 1988b ). Since these fields only 
exist for Cartesian and spherical geometry, the model of Title 
et al. (1993 ) could so far not be backed up by a self-consistent 
calculation. 



In a recent paper, Martens et al. (1990 presented a linear 
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force-free model for a fluted sunspot. The model is able to ex- 
plain some of the properties of fluted sunspots. Due to the high 
azimuthal wave number of the fluted force-free component, the 
scale height of this component is very small and the flutedness 
is confined to the chromosphere in this model. Furthermore, 
the magnetic loops corresponding to the dark filaments are very 
short in this model and the filaments do consist of a series of 
loops. This implies that the Evershed flow would have to con- 
sist of a phase-coordinated flow along many short loops. There 
is no explanation how this could be achieved. Therefore there 
is a clear need for improved models. 

The fundamental problem posed by the observations is that 
a genuinely three-dimensional solution of the magnetohydro- 
static equations is necessary. Though there has been some pro- 
gress concerning analytical solutions in three dimensions re- 
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cently (e.g. Neukirch 1995, 1997), these solutions do either not 
apply to the problem or they suffer from the same deficiencies 
as the linear force-free solutions. A different approach is there- 
fore necessary. 

In the present paper we present a method which allows us 
to find a self-consist ent ve rsion of the model originally put for- 
ward by Title et al. ( 1993 ). The method is based on an asymp- 
totic expansion of the three-dimensional force-free equations. 
The lowest order of this expansion has solutions which are for- 
mally equivalent to laminated force-free fields and a large num- 
ber of different solutions can be found. These solutions can be 
used to model the magnetic field of the penumbra of a fluted 
sunspot. 

In the following, we first briefly review the theory of lam- 
inated force-free fields and then present the expansion proce- 
dure. A specific solution which is useful as a model for a fluted 
penumbral magnetic field is presented and its properties are 
discussed especially in connection with a theoretical explana- 
tion of the Evershed flow. 



2. Basic Theory 

2.1. Non-existence of laminated force-free fields in cylindrical 
coordinates 

We assume that the magnetic field of the sunspot penumbra is 
force-free and therefore obeys the equations 



j x B = 0, 

V x B = noj, 
V • B = 0. 

From Eq. ([[]) we conclude that 

V x B = aB, 
B-Va = 0, 



(1) 
(2) 
(3) 



(4) 
(5) 



where we have absorbed into a. 

The observations (Degenhardt & Wiehr 1991 ; Title et al. 



1993; Rimmele 1995a) show that simple sunspots may have 



a fluted penumbra with rapid variations of the inclination of 
the magnetic field with azimuth (0), however with almost con- 
stant magnetic pressure at constant radius. Since the azimuthal 
component of the magnetic field is small compared to the other 
components and the sunspot field can be considered as being 
force free, also the electric currents seem to flow mostly in 
the meridional planes (0 ^constant). The appropriate model 
to describe such a magnetic field would be given by lami- 
nated non-linear force-free fields (Low 1988a, 1988b). How- 



ever, these fields do only exist in Cartesian and spherical coor- 
dinates; in cylindrical coordinates appropriate for the descrip- 
tion of sunspots such fields do not exist. This can be most easily 
seen as follows. The basic assumption is 



B = (B r ,0,B x ), 
a = a(4>). 



(6) 
(7) 



Then Eq. (g) is automatically fulfilled. If we take the curl of 
Eq. (^) and multiply the result by we get 



e • V x (V x B) = a 2 B 4 
On the other hand 

e f Vx(VxB) = -V 2 B 4 ~ 



Ida 

r d(f) 



x B = 0. 



2 0B r *2iBa 



2dB r 



r dcj) 



(8) 



■(9) 



It follows that B r may not depend on <f> which precludes fluted- 
ness for force-free fields in cylindrical coordinates. This means 
that only rotationally symmetric equilibria of this type are pos- 
sible and therefore it is impossible find an exact force-free so- 
lution appropriate for fluted sunspots. 

2.2. Asymptotic expansion method 

However, if we are mainly interested in the penumbral struc- 
ture, we could try to find an approximate solution with the 
desired properties, based on laminated equilibria. We start by 
writing Eqs. (^) and (Q) together with Eq. (||) in cylindrical co- 
ordinates: 



B, 



18_ 

r dr 



id_ 

r dr 
da 

7T + 
or 

(rB r . 



dBj, 
dz 
dB z 
dr 

T ) r 86 



I BB* 

r d(f> 
dB r 
dz 



1 da 

B Kd^ + B 

, IdB^ 



da 
dz 
dB z 



r dcj) 



aB r , 

aB<p, 

aB z , 

0. 

0. 



We now introduce a new radial coordinate 



2a 



(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



where a is a yet unspecified radius located somewhere inside 
the penumbra (actually a could be treated as a free parameter 
which can be used later to minimise the residual force resulting 
from the approximation). If the radial extent of the penumbra 
is 2R p , we now expand the equations above in the parameter 



Rpl a, 



(16) 



which we assume to be smaller than one (this is a kind of large- 
aspect-ratio expansion with 1/e being the aspect ratio). This 
type of expansion procedure together with the coordinate trans- 
formation Eq. (Q) has been introduced by Kiessling (1995) in 



the framework of toroidally confined plasma equilibria. 
With r = ay/1 + 2ew/R p we get 



1 



dB z 



dy/l + 2tzu/Rp 
dB r 
dz 



1 + 2ew/R 



= aB r , 
-dB z 



dm 



(17) 
(18) 
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dm 



'l + 2em/R p B cj 



1 



dB r 



ly/l + 2em/R p 



= aB z 



da 

1 + 2em/R p B r - — 
dm 



1 



da 



B<t> 7= = 7^=^77 + B z — — 0, 

ly/l + 2em/R p 9(p dz 



— ( x Jl + 2em/R pBr 



dB^ dB z 



0. 



ay/1 + 2em/R p dcj> dz 
We now introduce the following expansion scheme in e 

n— oo 

B r = Y, e n BW(w,4>,z), 

n=0 
n— oo 

n=l 
n— oo 

n=0 
n— oo 

a = J] e"a (,l) (n7,0,z). 



(19) 
(20) 

(21) 

(22) 
(23) 
(24) 
(25) 



n=0 



Note that is of order e whereas all other quantities are of 
order e° = 1. To lowest order in e we get: 



1 as 



(0) 



rf } dB z 0) 



dz dm 
1 dB r 0) 



/m day 3 ' /n\ 
Bi 0) — — + M 0) 



a d(f> 
da® 



0. 



(9c 



(92 



= 0, 



= 0. 



(26) 
(27) 
(28) 
(29) 
(30) 



dm dz 

Eqs. (E6h - (DO) are completely equivalent to the equations for 



laminated force-free fields in Cartesian geometry (Low 1988c) 
with the replacements x :— m and y : 

Multiplying Eq 
adding the two equations, we get: 

d 



by Bi 0) and Eq. @ by -B r ° } and 



/4 0)2 + ^ 0)2 ) 



o. 



(31) 



Eq. ( |29| ) can be solved by 

aW=a<°>(0). (32) 
Eq. (|30|) is solved by introducing a flux function A(i<7, 0, z): 



Inserting Eq. ( ]33[ ) into Eq. ( |27| ) we obtain 

Eqs. and (£H)) allow solutions of the form (Low |1988a i 

A(m, & z) = Re exp[i 7 (0)]) , (35) 

where £ = tu + iz. So in this approximation we have a field 
without a i?^ component that is potential in the nj-z-plane with 
an arbitrary modulation in the ^-direction. 

3. Modeling the penumbra of a fluted sunspot 

3.1. The Beckers-Schroter profile 

The solution (|5|) at first sight leaves a large degree of freedom 
to model fluted penumbrae. However, if one wants to calculate 
fields which reproduce the basic features of the observations 
the solution class is considerably restricted. We will illustrate 
this by considering a specific example. 

We assume that the basic data we have are the absolute 
value of the magnetic field (B(r)) and the inclination angle 
8(r,cj>,0) of the magnetic field with the vertical at the level 
z = 0. Within the framework of the theory developed so far, 
we can only represent spots which have little or no variation of 
field strength with azimuth. For reasons of mathematical con- 
venience, we assume that 



B(r) = 



Bo 



1 



{r/Lf 



(36) 



This is the radial dependence of the field strength deduced by 
Beckers & Schroter (1969) from their observations. Other ra- 



dial profiles like that of the Schatzman field (Schatzman 1965) 



could also be dealt with, but the Beckers-Schroter profile has 
the advantage of greater mathematical simplicty. 
From Eq. (|35|) the magnetic field is given by 



B r = Im 



B z = Re 



^ exp[ij((f))} 
— exp[i7(0)J 



The field amplitude is the given by 
B d ° 



d£ 



(37) 



(38) 



(39) 



With the coordinate transformation (15) we rewrite the 
Beckers-Schroter profile ( |3^ ) as 

Bo 



B(m) 



1 



2am 



(40) 



L 2 



Modulo a phase factor this profile corresponds to the absolute 
value of the complex function 



An 



4 



Neukirch & Martens: Self-consistent theory of fluted penumbrae 



on the real axis (z = 0). The most general complex function 
ip{£) with the same absolute value on the real axis can be writ- 
ten as 



m=FBs(o-Q(o, 

with 

\Q(a)\ = i 



(42) 



(43) 



on the real axis (z = 0). As we will see below the inclusion of 
an appropriate function Q is indispensible for a proper repre- 
sentation of the penumbral magnetic field. A simple example 
of such a function Q is given by 



Q(0 = exp(i<?£): 



(44) 



with q a real number. This example also illustrates that it is 
most convenient to write Q(£) in the form 



Q(0 = exp(i/(0). 



(45) 



The magnetic field components are the given by the expres- 
sions 



B,. 



BqL 2 exp(-Im(/)) 



2a 



L 2 



2a 



+ z 2 



L 2 



2a 



sin( 7 (0)+Re(/))- 

zcos( 7 (0)+Re(/)) 



(46) 



B, 



BqL 2 



exp(-Im(/)) 



2a 



2a 



L 2 



2 a 



cos( 7 (^)+Re(/)) + 

2sin( 7 (0)+Re(/)) 



(47) 



Another important piece of information given by the obser- 
vations is the inclination angle 5 of the magnetic field vector, 
i.e. the angle between the magnetic field vector and the vertical 
direction. For the magnetic field with the components given by 
Eqs. ( ^6| ) and <$t\), the inclination angle is given by the equa- 
tion 



5(r, cj), z) = 7 ((/>) + Re(/) — arctan 



, (48) 



\H7 + 



2a 



where the last term is the phase of Fbs- We remark that with 
the present theory, we can only represent fields that have an 
inclination angle which is additive in the dependence on the 
azimuth <f> and the radial coordinate w for 2 = 0. Therefore, 
we can write the inclination as 



6{tU, 4>, z) = 6 av {tu, z) + S mo d(4>), 



(49) 



with the identities 

6 av (w, z) = Re(/) — arctan 

Smod(0) = l{4>)- 



2 i „2 ' 



VO - 



L 2 +a 
2a 



(50) 
(51) 



The task now is to find a complex function /(£) which has the 
property Im(/) = for z = and which represents the ra- 
dial variation of the average inclination for z = 0. At first sight 
this might seem easy because we have already presented such 
a function above (/(£) = qO- A closer inspection of the in- 
clination angle and the magnetic field components for values 
of z reveals, however, that a naive choice of / will cause 
problems. For z = the phase of Fbs vanishes, but for z ^ 
it contributes to the inclination angle. Actually for z — > oo this 
contribution tends to —ir/2. Since the average inclination angle 
should be positive for z = this implies that the radial mag- 
netic field component will change its sign as z increases unless 
Re(/) also increases with z sufficiently fast. In other words, 
the field lines would bend backwards toward the sunspot center 
at a finite height z. This is of course not an acceptable solution. 
The real part of the linear function /(£) = q£ for example does 
not increase with z and therefore the linear function is not a 
viable choice. 

On the other hand, if Re(/) increases without bound for 
z — » oo the sign of both B z and B r will change an infinite 
number of times as the argument of sine, respectively cosine 
tends to infinity. This as well is not an acceptable solution. We 
therefore have to find a function /(£) which has a real part 
that tends to a constant as z goes to infinity and fulfills all the 
conditions mentioned previously. 

A function which fulfills all these conditions is 



/(£)=<? 



i-L a 



(52) 



where q, L a and Lb are real numbers. The real and imaginary 
parts of this function are 



Re(/) 
M/) 



= q 



{w - L a ){w - Lb) + z 2 
(m ~ L b ) 2 + z 2 
{L a - L b )z 
{■w - L b ) 2 + z 2 ' 



(53) 
(54) 



For z = the imaginary part vanishes as required and for vo 
fixed and z — > oo the real part has the limit q. Other func- 
tions may fulfill these requirements as well, but this form of / 
has the advantage of being relatively simple. Of course, if the 
radial variation of the inclination angle would be given by ob- 
servations, then the observations would prescribe the real part 
of / for z — 0. Here, however, we are merely interested to il- 
lustrate the method developed and therefore choose a simple 
form of /. Actually, once one has found one / with the desired 
properties, it is possible to add any other complex function / 
which is bounded for z — > oo and w fixed. An example of such 
a function is the linear function / = q£. This means that even 
with the restrictions discussed above one still has considerable 
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Table 1. Parameter values used in the example shown in Figs. 
0,|and| 



Parameter 


Meaning 


Value 


To 


outer penumbral radius 


1.0 


n 


inner penumbral radius 


0.5 


a 


origin of transformed radial coordinate 


0.8 


q 


amplitude of argument of phase factor 


7tt/9 




average inclination at inner boundary 


50° 




average inclination at outer boundary 


70° 


L a 


parameter of function / 


-0.562 


U 


parameter of function / 


-1.349 


&mod 


amplitude of variation 


18° 


III 


azimuthal wave number 


66 


B, 


field strength at inner boundary 


1500 G 


Bo 


field strength at outer boundary 


1000 G 


Bo 


amplitude of Beckers-Schroter profile 


1.8 


L 


length scale of Beckers-Schroter profile 


1.118 



freedom to model the variation of the average inclination angle 
with radius across the penumbra. 

3.2. An example 

In the present paper we restrict our treatment to the function / 
given in Eq. (52). We treat the parameter q as a free parameter 
which we choose such that we get a satisfying field line shape. 
The parameters L a and Li, are determined by imposing the (av- 
erage) inclination angles at the inner and outer boundary of the 
penumbra. This leads to two simple linear equations for L a and 

For the example we choose the outer boundary of the pen- 
umbra to be normalized (r — 1.0) and the inner boundary 
of the penumbra to be located at r; = 0.5. The origin of the 
transformed coordinate system is choosen as a = 0.8 locating 
it approximately halfway between the inner and outer boundary 
of the penumbra. The (average) inclination angles at the inner 
and outer boundary are chosen as S av .i — 50° and S av , = 70°. 
In this example we have chosen q = 7ir/9. The value of q 
influences the inclination of the field lines for z ^ 0. Higher 
values of q give a larger inclination for large z and vice versa. 
The resulting values for L a and L u are listed together with the 
other parameter values in Table |] 

For the modulation of the inclination angle with <j>, we 
choose a simple harmonic function 



5 mo d{4>) = l{4>) = 5 mo dS\n{m(, 



(55) 



The parameter 8 m od 1S the amplitude of the inclination variation 
which we choose to be 18° and m is the azimuthal wavenumber 
of the variation which we choose to be 66 in accordance with 



Martens et al. (1990- Of course, any other modulation with 
azimuth would be possible as well as long as it is periodic with 
period 2ir. 

The last two parameters to determine are those of the mag- 
netic field strength, namely Bq and L. These parameters will 
be determined by the measured values of the magnetic field 
strength at the inner and outer boundary (B. L and B a ). In the 




Fig. 3. Close-up of a part of Fig. ||. 



present example we have taken B t = 1500 G and B a = 1000 
G. Again, from the two conditions on the inner and outer 
penumbral boundary, we get two equations for the parameters 
Bo and L, which can be easily solved. For the field strengths 
used we also give the values of Bo and L in Table 111 The results 
are shown in Figs, [j] and |[ In Fig. [j] we show field line plots of 
the field lines in the planes <f> = 0, <fi = 7r/4m and <fi = 37r/4m 
corresponding to average, maximum and minimum inclination. 

In Fig. ^ we show a set of representative field lines having 
their foot points either close to the umbral-penumbral boundary 
or to the outer penumbral boundary. The field lines are chosen 
so that they outline locations of minimum and maximum incli- 
nation. The flutedness of the field is obvious. In Fig. |]we show 
a slice of the full field over two wavelengths in azimuth to em- 
phasize the details of the field structure. Especially in Fig. |3| 
we see that the field lines with the maximum inclination form 
long shallow loops across the penumbra as required for their 
association with the Evershed effect. 

We remark that the model could also be modified to i nclud e 
the return flux observed by Westendorp-Plaza et al. (1997). 



This could be done either by choosing different functions /(£) 
(or another set of parameters with the same function /(£)) or 
by extending the domain beyond the penumbral boundary and 
shifting the origin of the transformed radial coordinate system 
further out. 



4. Discussion 

4.1. Comparison with the linear force-free model 

There are clear differences between this model and the linear 



force-free model of Martens et al. (1996). In the present model 
the field strength drops off with height more or less like 1 / z 2 
for large z. The variation of the inclination angle with height is 
completely determined by the function / we choose, but in any 
case it will have a much weaker dependence on z than the linear 
force-free model. In the linear force-free model, the magnetic 
field drops off exponentially and the scale height of the fluted 
part is very small. This is actually one of the shortcomings of 
that model. 

Another difference is that the parts of the penumbra where 
the field is nearly horizontal do not form rather short loops but 
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0.50 0.75 1.00 0.50 0.75 1.00 0.50 0.75 1.00 



Fig. 1. Field line plots in the planes </> = corresponding to average inclination (left), <f> = ir/4m corresponding to maximum 
inclination (middle) and <f> — Sir /Am corresponding to minimum inclination (right). The field lines plotted cross the lower or 
left boundaries at the same location in all three plots to facilitate the comparison. Therefore the distance between field lines does 
not necessarily reflect the strength of the magnetic field. Though not shown here, we emphasize that it is possible to extend the 
model in such a way that the return of the low-lying field lines to the surface is included in accordance with the observations by 
Westendorp-Plaza et al. (|1997|). 



quite extended radial loops which can more easily account for 
the observed Evershed flow. 

Whereas in the linear force-free model the direction and the 
amplitude of the current density are determined by the mag- 
netic field because a is a constant, in the present model is 
a function of the azimuth. The relation between and 7 can 
be calculated either from Eq. (Efl) or Eq. (E8). We get 



1 d 7 

a d(f> 



With the form of 7 given in Eq. (|55j), we obtain 



-8 m od cos(ra</>) 



(56) 



(57) 



Since the cosine has its extrema where the sine has its zeros 
and vice versa, we have the situation that the maximum current 
density flows along field lines having the average inclination, 
whereas the current density vanishes along field lines having 
the maximum or minimum inclination. Furthermore the direc- 
tion of the current flow changes its sign every half wave length 
in azimuth. T his co incides exactly with the discussion given 
in Title et al. ( 1993 ) and sketched in their Fig. 17. This vindi- 
cates our approach because it was our aim to come up with a 
self-consistent version of their schematic model. 



4.2. Quality of the approximation 

An important point to investigate is the quality of the approxi- 
mation scheme that we have presented. A good way to do this 
is to investigate the magnitude of the residual force due to the 
approximation. To be able to judge the quality of the approxi- 
mation, we need to compare the residual force to a quantity of 
the same dimension. A convenient measure for the strength of 
the force is B 2 (r, z)/hqL. We then obtain 



\j x B\ 
\B\ 2 







r --l 
a 


dB z 

dm 


\B\ 




\B 





(58) 



Note that the residual force has only a poloidal component. One 
can immediately see that the residual force vanishes at the ori- 
gin of the transformed coordinate system (r = a). In Fig. |] 
we show a surface plot of the residual force as function of r 
and z for = 0. The plots for the locations of minimum and 
maximum inclination do not differ very much from this plot. 
It can be seen that the residual force is small close to the ori- 
gin of the transformed radial coordinate system and then rises 
almost linearly in r away from that origin. This can be under- 
stood by an inspection of Eq. (|58[). Most of the variation of 
the residual force expression obviously comes from the factor 
\r/a — 1|. The factor dB z / dzu / B just seems to be a minor 
modulation of the first factor. The plot also shows the limits of 
the approximation scheme. At the boundaries (in r) the dimen- 
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Fig. 2. Three-dimensional plot of sets of field lines originating close to the inner and outer penumbral boundary. The flutedness 
of the field is obvious. The azimuthal wave number is m — 66 in this case. 




Fig. 4. Surface plot of the relative strength of the residual j x 
.B-force above the r-z-plane for = 0. The magnitude of the 
residual force stays well below 1, indicating that the quality of 
the approximation is reasonable. 



sionless residual force has reached a value of almost 0.6 which 
is at the very limit of what can be considered as acceptable for 
a small parameter. On the other hand, we have only considered 
the lowest order of the expansion scheme and higher order cor- 
rections could make the representation of the field even better. 
One should also keep in mind that by an expansion procedure 
like this one can usually not expect to get a convergent series 
but only an asymptotic series (though in a mathematically rig- 
orous sense we would still need to prove that the series is in- 
deed asymptotic). 

We mention as a possibility that in principle, one could inte- 
grate the amplitude of the force over the volume under consid- 
eration and minimise this integral with respect to the parameter 
a. This would give a kind of optimum value for the location 



of the origin of the transformed r-coordinate. Since we expect 
that the increase in accuracy achieved by such a procedure will 
be small we have not carried this out here. 



5. Conclusions 

We have presented an expansion scheme that allows us to cal- 
culate self-consistent force-free solutions of the magnetohy- 
drostatic equations with the basic properties of fluted sunspot 
penumbrae. The lowest order of the expansion produces equa- 
tions which are mathematically equivalent to the equations 
describing laminated force-free solutions in Cartesian coordi- 
nates. The magnetic field then has no toroidal component in 
the lowest order and the magnetic field strength does not vary 
with azimuth. The field inclination, however, can have an arbi- 
trary variation with azi muth. I n comparison with the force-free 
model by Martens et al. (1996 ) the scale height of this variation 
is generally much larger. 

We have calculated explicit solutions which have a pro- 
file of the magnetic field strength in the plane z = t hat ex - 
actly matches the profile given by Beckers & Schroter ( 1969 ). 
The field lines emerging from the plane z — extend across 
the penumbra in the radial direction and for the maximum 
inclination angle stay relatively shallow. Our model provides 
for long, low-lying radial loops that return to the photosphere 
just beyond the boundary of the penumbra, as is required for 
the Evershed flow, and confirmed by recent observations of 



Westendopr-Plaza et al. (1997) 
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